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ABSTRACT 

Preliminary results of an experiment to determine the 
accuracy and the effects of varying amounts of environmental 
data at the NAVTORPSTA, Keyport, Washington, Dabob Bay 
Range are presented. A 75-KHz source was tracked at six 
depths at each of six horizontal ranges. The data were 
analyzed by the NAVTORPSTA NUTRACK III isovelocity-layer 
computer program and by an isogradient-layer program 
developed here, STUTRACK I. The results show that a low 
signal-to-noise ratio is not a problem even at long distances 
and that an isogradient analysis yields superior results 
with less environmental data than does the isovelocity 
method. 
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I . INTRODUCTION 



A. DABOB BAY FACILITY 

A facility for the three-dimensional tracking of 
surface ships, submarines, and torpedoes is located at 
the Dabob Bay facility of the Naval Torpedo Station 
(NAVTORPSTA) , Keyport, Washington. Chosen because of its 
favorable oceanographic features and its proximity to the 
NAVTORPSTA, Keyport, Dabob Bay has been used for torpedo 
testing since 19^9. It consists of a 250-KHz "high 
frequency" installation no longer operating, and a 75-KHz 
"low frequency" primary tracking system enveloping nearly 
the entire bay providing for underwater tracking in a 
volume approximately 30,000 feet by *1,500 feet by 600 feet 
deep . 

Dabob Bay, a branch of Hood Canal, is in a deep (600 feet) 
depression adjacent to the Olympic Mountains. The Quilcene 
River flows into the northernmost part and the Dosewallips 
River into the southernmost. The sides of the bay are 
precipitous and predominately rocky, while the bottom is 
mud. Annual precipitation averages about 51 inches and 
while average snowfall is about 15 inches, any one snowfall 
seldom remains longer than a few days. Winds generally blow 
along the length of Dabob Bay; i.e.. North and South. 

Southerly winds are usually accompanied by clear skies and 
occur most often in summer. Normal wind velocities are 
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5 to 15 miles per hour. Tidal levels in Dabob Bay range 
from about -5 to +15 feet. 

Water temperatures change markedly from season to season 
and even from day to day. The changes are most noticeable 
near the surface but may be observed even at greater depths. 
Salinity varies considerably with depth as fresher surface 
water from terrestrial runoff overlays more saline water. 

The variation is also quite seasonal, reaching a maximum in 
the spring. The main contributor to sound velocity changes 
is water temperature, and the resulting sound velocity 
profile is representative of the temperature gradient. 

B. MOTIVATION 

NAVTORPSTA, Keyport, has developed a computer program, 
NUTRACK III, which divides the bay into laterally homogeneous 
isovelocity layers and traces a sound ray from a three- 
dimensional hydrophone array mounted just above the bottom, 
back to its source. 

An experiment was undertaken to provide comparison of 
the three-dimensional position of a source with respect to 
a receiver as determined by various ray-tracing techniques. 
Prom this comparison it was hoped that an optimum ray-trace 
procedure would be indicated and that a better understanding 
of the significance of environmental variations would result. 
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II . PROCEDURE 



A. DETAILED RANGE DESCRIPTION 

The basic underwater tracking system on the 75-KHz 
range is comprised of three components: hydrophone array, 

computer system, and vehicular instrumentation. The ship- 
board transducer emits a 75-KHz acoustic pulse in synchronism 
with a master clock at the computer site. This pulse is 

detected by each of the four hydrophones (R ,R ,R and R ) 

x y z c 

in the array and mixed with a local oscillator having a 
frequency unique to that hydrophone. The outputs of these 
oscillators are preamplified and fed to a multiplexer, which 
sums the signals and transmits the composite signal via an 
underwater cable to the computer site, where it is processed 
to determine the in-water transit time to each hydrophone 
<T x> T y) T z , and T 0 ). 

1 . Hydrophone Array Configuration 

The hydrophones in the NAVTORPSTA Range are arranged 
in a short-baseline system in groups of four, each one 
located on each of four adjacent corners in an imaginary 
30 -foot cube, thereby defining the orthogonal coordinate 
system in which the measurements were made. These hydrophone 
arrays used buoyant spheres which exert an upward force in 
excess of two thousand pounds to keep them upright. The 
array is not free to rotate about the Z axis, but may tilt 
as much as two degrees from the horizontal in strong currents. 
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This tilt is measured by servo-pendulum transducers whose 
specified resolution is 0.00833 degrees, accuracy 0.025 
degrees, and linearity 0.05 percent over the full scale. 

The Range consists of six arrays (Pig. 1) aligned 
along a line 450 feet east and parallel to the Range axis, 
spaced at 6000-foot intervals, with the 00 array being the 
farthest north. The Range axis is 191° 1 8 * 14.1" T. 

2 . Computer System 

The computer system, a Scientific Data Systems (SDS) 
model 920, is installed on Zelatched Point and consists of 
three main subsystems: the signal processing subsystem, 

the data collection subsystem, and the computer subsystem. 

The signal processing subsystem, the link between the hydro- 
phones and the data collection subsystem, processes the 
multiplexed signal received from the arrays, discriminates 
against unwanted signals, and determines the in-water transit 
time to each hydrophone. The data collection subsystem 
establishes the master clock timing by which all timing is 
computed. It also prepares all the array signals for the 
computer subsystem, which calculates the tracked object’s 
position, prints the tracking data, and records the data on 
magnetic tape. For real-time examination, the data are 
plotted on X-Y plotters. 

B. EXPERIMENTAL PROCEDURE 

A 92.5 dB acoustic source with approximately ten pounds 
of ballast was suspended from the stern of the sound boat 
to 25, 50, 75, 100, 150, and 200-foot depths at six 
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horizontal ranges (stations) from 900 feet to 4700 feet 
from the 02 hydrophone (Fig. 2). The first station was 
located approximately 900 feet north of hydrophone array 02, 
and 450 feet east of the Range axis. Successive stations 
were located parallel to the axis of the Range; the second 
being directly over the array, the third 600 feet past it 
(diametrically opposite the first station), and then at about 
1500-foot intervals to about 4500 feet, beyond which signal- 
to-noise problems made further examination impossible. The 
source emitted one 1.3 millesecond 75-KHz pulse per second 
for sixty to ninety seconds at each depth at each station. 

The position of the boat on the Range was found in the 
horizontal from the Range’s Autotape system, a portable and 
commercial microwave system used exclusively for ship 
tracking that yields accuracies of ±3 feet. The position 
of the acoustic source with respect to the Autotape antenna 
was determined from the length of cable let out, the 
horizontal distance from the antenna to the stern of the boat, 
and the boat’s heading. Tidal data were recorded and used 
for correcting depth coordinates to a mean low-low water 
(0-foot tide) reference plane. The array is known to be 
585 feet below MLLW, so that the depth of the source below 
the surface and the tide measurement allows calculation of 
the vertical distance between the source and array. All 
measurements were made with respect to the array's coordinate 
system, although the array's orientation and location on the 
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Range were known, making it possible to transform this data 
to a true position on the Range. The X axis of the 02 array 
was 25° 24' 36.0" T. 

A NAVTORPSTA digital velocimeter that has an accuracy to 
±1.0 ft/sec and resolution of ±0.7 ft/sec measured the sound- 
velocity profile at each station. NAVTORPSTA instrumentation 
at Dabob Bay measured and recorded the acoustic travel times 
from the source to each of the four hydrophones in the array 
to six decimal places for every source depth at each station. 
The three-dimensional position of the source with respect to 
the array was calculated from these travel times by the 
NAVTORPSTA NUTRACK III computer program. 

The data were gathered on 21 March, 1974 during a neap 
tidal period and while the sea surface was smooth and glassy, 
although there had been ripples earlier in the morning. The 
day was bright and sunny with scattered, thin, high overcast, 
the temperature was 55°F 3 and no noticeable wind was present. 
The shiny stainless steel transducer disappeared from view 
when about eight feet deep, indicative of thick biological 
concentrations (gelbstaffen) in the water. Rudimentary drift 
measurements of boat movement relative to surface water were 
taken by throwing paper into the water and watching its 
relative motion with respect to the boat, which indicated 
a value of approximately 0.3 ft/sec. 
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III. THEORETICAL DEVELOPMENT 



The computer program analysis includes corrections for 
the fact that acoustic rays are not straight lines in a 
non-homogeneous medium. Two simple ray-tracing techniques 
exist: the water mass divided into layers of isovelocity 

water or into layers with a constant velocity gradient. 

The isovelocity method is employed in the NAVTORPSTA 
NUTRACK III program; while employment of the isogradient 
technique is one of the purposes of this present experiment. 

A. NUTRACK III 

1. Raw Array Coordinates 

The geometry involved in the calculations of the 
raw array coordinates is shown in Fig. 3* As described 
earlier, the four hydrophones are on four adjacent vertices 
separated by a distance d, along the edge of the cube. The 
origin, 0, of the array coordinates is at the center of the 
cube with the orthogonal coordinates parallel to its edge. 
Let the transit times of the acoustic pulse from the 
tracked source to the c hydrophone and the Y hydrophone 
be T and T , respectively, then 

(cT c ) 2 = (RAWY + |) 2 + (RAWX + |) 2 + (RAWZ + |) 2 (l) 

and 

(cT y ) 2 = (RAWY - |) 2 + (RAW X + |) 2 + (RAWZ + |) 2 . ( 2 ) 
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2 2 

Taking the difference between T and T , expanding, 

c y 

and combining terms yields 



RAWY = (T c +T y )(T c -T y ) . 



(3) 



In similar fashion, analogous expressions may be 
obtained for the X and Z coordinates : 



RAWX = ~ (T +T )(T -T ) 
2d c x c x 



(4) 



and 



MWZ = (T 0+ T z )(T c -T z ) 



(5) 



where 

c = the average sound velcoity, about 4860 ft/sec; 
d = the distance between hydrophones, 30 ft; 

T c ,T x ,Ty,T z = travel times in seconds. 

The travel times are determined from a clock whose 
period is 0.4 x 10 ^ sec. (2.5 mHz), which is sufficiently 
high to measure precisely the arrival time of a 75-kHz 
signal (about 8 counts per quarter period). A correction 
factor must be applied to correct for the array receiver 
integration (pulse width descrimination) . Let 
I = the integration time, 800 ys 

T^ = the measured travel time to the R c hydrophone; 

T = actual travel time, 
c 

so that 
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i 



and 



T ’ = T + I 

c c 

t ; ■ T y + 1 > 

T c + T y = (T i + T y ' 2I) > 

T c - T y “ T i * ' 



The relations for the raw coordinates are now 



RAWY = fj (T’ +t^-k)(t;-t^) , 



( 6 ) 



where K = 21 
t ime s . 



1.6 x 10 sec and T^ and T^ are the measured 



Similarly 



RAWX = §-r (T* + T' - K) (T ’ - T' ) 
2d c y c x 



( 7 ) 



and 



RAWZ = — (T ' + T ' - K ) (T 1 -T') 
2d c z c z 



( 8 ) 



2 . Tilt Correction to Determine Corrected 
Array Coordinates 

The raw coordinates previously computed are 
referenced to a non-horizontal reference plane. The array 
is instrumented to sense the angles the array axes make 
with the horizontal, c&lled XTILT, or rotation about the Y 
axis, and YTILT, or rotation about the X axis, respectively. 
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These are given from 



XTILT = (TILTX - 505) x .00833 degrees 

and 

YTILT = (TILTY - 505) x .00833 degrees 

where TILTX and TILTY are sensor outputs and 505 is an 
arbitrarily determined sensor reading for a perfectly 
horizontal orientation. 

The correction of RAW coordinates for array tilt 
is merely a matter of geometry. Figure 4 shows that 

CORX = [RAWX] [COS (XTILT) ] = [RAWZ ] [SIN (XTILT) ] . (9) 

Since XTILT is small, it may be assumed that COS(XTILT) = 1. 
The relationship then becomes 

CORX = RAWX - [RAWZ ] [SIN ( XTILT ) ] . (10) 

Similarly 

CORY = RAWY - [RAWZ ] [SIN (YTILT) ] (11) 

and 

CORZ = RAWZ + [RAWX] [SIN (XTILT) ] + [RAWY ] [SIN (YTILT ) ] 

(12) 

These relationships are valid for both positive and negative 
tilt as the SIN term also changes sign to maintain the 
relationship . 
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3. Ray Path Refraction 

In a homogeneous body or water, the minimum time 
path between an output transducer and receiving hydrophone 
is a straight line. When the speed of propagation is 
spatially varying due to temperature and salinity variations, 
refraction occurs and the minimum time path is no longer a 
straight line. Corrections must then be made to determine 
the true position of the source. 

Prom the classical laws of Fermat and Snell, a 
minimum-time path is computed, working back in time from the 
array center in the direction of the "apparent" position and 
tracing the array through each velocity layer. The process 
is continued upward in space and backward in time until the 
measured travel time is consumed. To determine ray-path 
refraction, the sound speed is taken at one-foot increments 
and averaged over ten feet. Since the array depth is 
585 feet, there are 58 ten-foot layers and one five-foot 
layer. The vertical angle of entry of the ray at the 
array, computed from the tilt-corrected array coordinates 
(Fig. 5) j is 



A = sin- 1 r cop . ^ ] 1/2 . (13) 

° CORX 2 + CORY^ + CORZ 2 

The entry angle A o computed from the CORR array 
coordinates is assumed to be equal to the entry angle for 
the minimum time path. This assumption is good if the 
source is far enough so that the wave fronts are plane over 



20 



the dimensions of the array, and the speed of sound is 
constant over the dimensions of the array. 

The travel time, T, to the origin is computed from 



ir m ir > (l4) 

c c 

where 

S = [CORX 2 + CORY 2 + CORZ 2 ] 1/2 (15) 

and 

R„ = [(CORX + 15) 2 + (CORY + 15) 2 + (CORZ + 15) 2 ] 1/2 , 

V-r 

(16) 



where R is the slant range to the C array hydrophone, 
c 

and S is the slant range to the array center. 

It may be seen that the assumption of an average 
velocity, c, in computing RAWX, RAWY, RAWZ is not sensitive 
as the ratio of S to R q is what is really important and any 
small error will tend to cancel out. 

The horizontal distance is 

H = [CORX 2 + CORY 2 ] 1/2 . (17) 

The first layer processed is between the deepest 
depth in the velocity profile that is shallower than the 
array, and the array itself; 

AZ 1 = A d " Z 1 5 (l8) 

where A^ is the array depth. 
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The incremental slant range is 



AZ X 

AS 1 = SIN 

o 



( 19 ) 



The incremental travel time in this layer, used to increment 
the pulse travel time T, is 




( 20 ) 



where c-^ is the average velocity in the layer. This 
incremental travel time is used to decrement the pulse 
travel time T. The horizontal distance for this time is 



AH, = AS, COS A . (21) 

1 1 o 

The values AH and AZ are accumulated as the time clock is 
decremented layer by layer. 

Using Snell's law, the change in angle due to 
refraction upon entry into a new velocity layer is 

COS A c, 

o _ 1 

C0S(A q + AA) c 2 



or 



COS A 1 



COS A — 
° C 1 



( 22 ) 
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where 



A 1 = A o + AA 

and ' 

SIN A ± = (1 - COS 2 A 1 ) 1/2 
The thickness of the layer is 



( 23 ) 



AZ 



2 




( 24 ) 



This process is repeated and the original travel 
time T is decremented with the incremental layer travel 
times until T is consumed. The ray-bent X, Y, and Z 
( RBNX , RBNY, RBNZ ) are 



RBNX = ( 



CORX . 
H } 



I AH 



( 25 ) 



RBNY = ( ) £ AH 



( 26 ) 



and 



RBNZ = £ AZ 



( 27 ) 



This raypath migration process assumes a laterally 
homogeneous medium within the neighborhood of each array. 
Independent velocity profiles are not taken at each array. 
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This will introduce small discontinuities in track as the 
tracked vehicle is passed from one array to another. 

4 . Error Estimation 

Returning to Eq. and assuming equal errors, 

AT = AT = At in time measurement; the maximum expected 

C X 

error in RBNX (AX) is found from 



RAWX + AX = §-r [(T +At) 2 -(T -At) 2 ] 

c. Q C X 




C T c- T x] + %- CT e At+T x At] 



( 28 ) 



to be 



AX - [T c + T x ]At . 



At large ranges, T « T , so 

C X 



= IsL 

d J ' c L 



AX = T At , 

r* * 



but c*T = R , so 
c 



¥ * 2 <§> At • 



Similarly , 



AY i 2 (— ) At 
R 



and 



R " 2 (f) At • 



( 29 ) 



( 30 ) 



( 31 ) 

( 32 ) 



2 'I 



Assuming the maximum error in time is one-fourth of the 

period of the 75-kHz signal. At a 3-3xl0 - ^ sec, so 

c 



AX 

R 



AY 

R 



AZ 

R 



* 1 x 10 
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B. ISOGRADIENT PROGRAM "STUTRACK I" 

An isogradient computer program, STUTRACK I, was 
developed for comparison with NUTRACK III. 

1 . Linear Gradient Theory 

Figure 6 shows a highly exaggerated speed of sound 
profile which has been approximated by two layers each 
possessing a constant gradient. In each layer, the speed 
of sound may be given as 



c ( Z ) = c Q + GAZ (33) 

•s. 

where G is the gradient; G is negative in the upper layer 
and positive in the lower layer. Applying Snell's Law for A, 



- C - ° S - - = CONSTANT , ( 3 h ) 

c 

one may see that the raypath curves downward in the upper 
layer and upward in the lower layer. Let the radius of 
curvature for the arc between the angles A^ and A^ in the 
upper layer be R. Then 
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( 35 ) 



AZ = R(COS A, - COS A~ ) = R COS A- ( — -) . 

1 d. . J C ^ 

However, from Eq. 33, 

— c 2 = — GAZ 
and with Eq. 34, 



R = 



zi e- 
G COS A 



(36) 



Thus for each raypath the value of R is a constant 
in a linear gradient layer, and the path is an arc of a 
circle. The center of the circle lies at the depth for 
which A would become 90° if the sound profile for the upper 
layer were to descend to greater depths. 

Assuming that the initial angle of depression of the 
raypath is A Q , from the geometry of Fig. 6 and with Eq. 36, 
the change in depth and horizontal range become 

AZ * - G COS A (C0S A o - C0S A> (37) 

o 

and 

Ar = - g~c5§' T ~' (SIN A - SIN V ’ (38) 

where Z is the increase in depth and r is the increase in 
range from the reference position, and A is the angle of 
depression of the raypath at this range and depth. 



26 



When AZ is known, the new angle A may be found by 



A = COS 1 [COS A (— + 1)] 

o c o 



( 39 ) 



The time-of-f light along the raypath may be 
simplified to 



At 



A 

-1 , dA 

G A J COS A 



(40) 



which upon integration yields 



A j_ 1 COS A ^ / 

At - g ln[( 1 + SIN A ) ( 



1 + 



SIN A n 



COS A 



)] 



(41) 



2 . STUTRACK I Construction 

STUTRACK I (App . A) is designed to allow the deter- 
mination of a minimum amount of environmental information 
necessary for agreement with the NUTRACK III results. The 
program determines the tilt-corrected coordinates and the 
raypath entry angle into the array in the same manner as 
NUTRACK III, but uses Eq. 39 to calculate the exit angle A 
from the layer, Eq. 40 to calculate the horizontal range 
traversed, and the time-of-f light from Eq. 4l. The thickness 
of the layer, AZ, is determined on a point-to-point basis 
from the depths from which the velocity profile was established. 
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The values Ar and AZ are accumulated and the original travel 
time T is decremented as in NUTRACK III until 

TG > T + e , 



where 



and 



TG = E At 



e - . 0005 sec . 



This value of e ensures an accuracy of E AZ and E Ar 
within ±3 inches. If TG is within e of T, the process 
stops and the X, Y, and Z coordinates for this isogradient 
technique (XG,YG,ZG) are calculated as in NUTRACK III, 



and 



XG 

YG 



CORX 

H 



CORY 

H 



E Ar , 
E Ar , 



ZG 



CORZ 

H 



E Ar 



(42) 

(43) 



(44) 



When 



TG > T + e , 

the program returns to the previous layer and a bracket- 
and-halving process begins in AZ until the proper AZ is 
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found to satisfy the e criterion in time. XG, YG, and ZG 
are then found by Eqs. 42, -# 3 , 

STUTRACK I is also designed to permit the utilization 
of the velocity-profile data in varying incremental step 
sizes to allow a determination of the effect of a varying 
amount of environmental data. For example, if velocimeter 
data is gathered in two-foot steps, an incremental step size 
(INC) of 1 would automatically set AZ = 2 ft, and an INC of 
10 would set AZ = 20 ft. Consideration of the source depth 
must be made prior to the invocation of STUTRACK I in its 
present form, as the setting of too large an INC can create 
E AZ greater than the array depth of 585 feet, resulting in 
an error message. 

STUTRACK I is also designed to analyze the data with 
the same isovelocity technique as NUTRACK III, but in 
varying incremental step sizes for comparison with the 
isogradient method. These output coordinates are XI, YI, 
and ZI. 

The inputs into STUTRACK I are: velocity profile, 

increment size, tilt correction, time of day in hours, 

minutes and seconds, T , T , T , T , RBNX, RBNY, and RBNZ. 

x y z c 

The output is six histograms of the differences between the 
NUTRACK III, ray bent coordinates, and their respective 
values as determined by STUTRACK I using both the isogradient 
and isovelocity methods. These differences and the STUTRACK I 
values are also printed. 
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IV. EXPERIMENTAL RESULTS 



A. VELOCITY PROFILE 

Figures 7 and 8 Illustrate that the velocity profiles 
for station 1 (taken at 1115 hrs) and station 6 (taken at 
1315 hrs) are equivalent below 85 feet. (See Fig. 2 for 
station positions.) The data for station 6, represented by 
the "+"'s, shows the "afternoon effect" of solar heating. 
Both profiles increase in velocity with depth at a greater 
rate than would be caused by pressure alone, implying a 
similarly increasing temperature profile. The colder 
surface water was fresher, however, due to the heavy spring 
rainfall and runoff from the snow-capped mountain peaks. 

B. AUTOTAPE DATA 

The Range's microwave Autotape system was used to track 
the sound boat. Its precision was much less than expected 
so that comparison with Range data was not attempted at this 
time and will have to wait until the Autotape data is 
subjected to extensive analysis. 

C. ANALYSIS OF NUTRACK III RESULTS 

1. One-Way Transmission Loss 

The one-way transmission loss (TL,) is determined by 
an assumption of spherical spreading and the relationship 
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TL(dB) = 20 LOG R + aR 



> 



( 45 ) 



where 



a = attenuation constant in sea water 
for 75-kHz sound. 



- 0.03 dB/meter. 



Figure 9 shows that at short ranges (station 1) the TL is 
dominated by the spherical spreading term. Beyond a "critical 
range" of approximately 1000 feet for 75-kHz sound, the 
losses due to absorption of sound in sea water will overpower 
the spreading and the transmission loss will increase 
rapidly at greater distances. The TL from station 6 to the 
hydrophone array is approximately 105 dB, but is only 57 dB 
from station 1 . 

2. RBNX, RBNY , RBNZ and Slant Ranges 

For Stations 1-25 and 6-50 

Figures 10-13 show the 80 data for station 1-25 
plotted versus time. Figures 14-16 represent the same data 
as Figs. 10-12 but with a linear drift correction. Figures 
17-21 are similar curves for the 66 data for station 6 - 50 , 
except that a drift correction was applied only to the slant 
range curve. These two stations were chosen for preliminary 
investigation since station 1 was closest to the array while 
station 6 was farthest away. The RBNX, RBNY, and Slant Range 
scales for each station are nearly equivalent, but RBNZ data 
are plotted in a 5/1 ratio, since station 6 was five times 
the slant range of station 1. In all cases, station 1 data 
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show good point-to-point continuity while station 6 data, 
except for Slant Range, are very erratic. Microstructure 
was examined as the cause of this variation using the 
procedure of Medwin [2], but was determined to be capable 
of producing a maximum fluctuation of only ±1 foot, or only 
a small fraction of the actual variation. 

3. Transformation to Spherical Coordinates ' 

A transformation from the orthogonal coordinates 
RBNX, RBNY, and RBNZ to a spherical system was achieved 
using the following relations, 

R = (RBNX 2 + RBNY 2 + RBNZ 2 ) 1/2 , (46) 



$ 



SIN 



( RMZ) 






and 



0 



TAN 



- 1 / RBNY v 
' RBNX ' 



5 



where $ is the angle of elevation looking from the array to 
the source, and 0 is an arbitrarily defined azimuthal angle. 
A graph of $ versus 0 is a scatter diagram of the projection 
of the data onto the surface of a sphere of radius R. 

Figures 22-29 are station 1-25 plots of $ and 0 as 
calculated from: 

a. Ray bent coordinates (Figs. 22-24), 

b. Ray bent coordinates, amplified, and with linear 
drift correction (Figs. 25-27), 

c. COR coordinates, amplified (Figs. 28 and 29). 
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Figures 22-24 show a consistent, nearly linear drift, with 
good point-to-point continuity. Figures 26 and 29 are 
identical, as should be the case since NUTRACK III assumed 
a laterally homogeneous medium, thereby keeping the ratio 
Y/X constant. Figures 25 and 28 show similar but not 
identical variations, which was also expected since the COR 
coordinates assumed a vertically homogeneous medium and the 
raybent coordinates consider the inhomogeneities. 

Figures 30-32 are station 6-50 plots of $ vs. time, 

0 vs. time, and $ vs . 0, respectively, as calculated directly 
from the raybent coordinates, and plotted to the same scales 
as Figs. 25-27. The variations are approximately the same 
in magnitude as for the same plots at station 1-25. However 
the scatter diagram for station 1-25 transforms into a five- 
foot square patch on the surface of a sphere of radius 900 
feet. At station 6-50 this uncertainty grows with the ratio 
of the radii to a patch twenty-five feet square. 

4 . Error Determination 

An error determination was made from the data from 
Figs. 12, 14, 15, and 17-19 by subtracting each datum from 
the average value drift-corrected coordinate and dividing 
by R. Variations from the averages were found to be less 
than 5x10 with a standard deviation (a) less than 10 . 

This latter is comparable with the theoretical error deter- 
mination of section III-A-4 on the basis of an error in 
determining the arrival time which is independent of range. 
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D. COMPARISON OF NUTRACK III TO STUTRACK I 



Comparison of STUTRACK I data in both the isogradient 
and isovelocity modes with the NUTRACK III raybent coordinates 
was achieved by the following relationships 

DXG = RBNX - XG (1*7) 

DYG = RBNY - YG 
DZG = RBNZ - ZG 

and 

DXI = RBNX - XI (48) 

DYI = RBNY - YI 
DZI = RBNZ - ZI . 

Station 6-50 data (Figs. 36 - 38 ) are consistent with 
station 1-25 data (Figs. 33-35) s but much more dramatic 
since it was taken at a much greater range. These graphs 
are plots of DX, DY, and DZ for both the isogradient and 
isovelocity cases as a function of incremental step size. 

They show that the STUTRACK I results are essentially 
equivalent to the NUTRACK III data when INC = 1. However 
as INC is increased, the isovelocity technique gives 
consistently poorer results while the isogradient method's 
results remain adequate as INC increases. In all cases the 
STUTRACK I coordinates followed the NUTRACK III raybent 
coordinates point-for-point , but were biased to one side. 

A few seconds of computer time was also saved as a result 
of the isogradient technique at a greater INC than for the 
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comparable calculations via the isovelocity method since 
fewer layers had to be examined to ’achieve the same, or 
better, results. 

Figure 39 is a histogram of DXG, INC = 1, at station 
6-50. Its Gaussian nature is characteristic of all the 
histograms of DX, DY, and DZ at every INC tested, at both 
stations . 
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V. CONCLUSIONS 



A. NUTRACK III DATA 

The large variations in RBNX, RBNY, and RBNZ compared to 
the smoothness of the Slant Range may be explained by an 
examination of Eqs. 3, 4, 5, and 15. The individual 
coordinates are dependent upon the difference of T and T , 

C X 

T, or T . Since this difference is small, a small error 
y z 

in time will make a large fractional error in the computed 
coordinate. The uncertainty in signal discrimination, 
addressed in section III-A-4 and estimated as 3-3 x 10~^ sec, 
becomes important, causing significant random errors from 
pulse-to-pulse . This error becomes very large for any 
position for which the wave arrives at the two hydrophones 
at nearly the same time. For calculation of the Slant Range 
the arrival time is used alone so the fractional error is 
much less. 

B. COMPARISON OF STATION 1-25 TO 6-50 

Comparison of station 1-25 data to that from 6-50 shows 
the latter to have greater spatial variability. The trans- 
formation to spherical coordinates shows this difference is 
not present in the angular display. This is what would be 
expected if the time errors were the important sources of 
variability. This conclusion is upheld by an experimental 
error determination at both stations that agrees with the 
theoretical error induced by the one-fourth period 
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uncertainty in the signal arrival times. Out to the 
greatest distances studied (over twice that used in normal 
Range operations), the precision in measuring arrival times 
is unchanged even though the signal-to-noise ratio decreases 
by nearly 50 dB. Therefore, it is concluded that for 
conditions similar to those prevailing at the time of this 
experiment any attempt to improve the signal-to-noise level 
by increasing the source level will be nonproductive and 
wasteful. Improved Range accuracy can only be obtained by 
improving the accuracy of the time measurements (above its 
present value of about one-quarter period of the 75-kHz 
signal) or increasing the baseline of the hydrophone array. 

C. EFFECT OF ENVIRONMENTAL DATA 

The graphs of DX, DY, and DZ for both stations show 
conclusively that for the conditions prevailing at the time 
of this experiment the isogradient technique of STUTRACK I 
provides better results with less environmental data than 
does NUTRACK III. This is consistent in each of the plots 
but. is explicit in DZ at station 6, where the isovelocity 
layer method results in a linear growth in DZ as the 
incremental step size (INC) is increased, but little change 
in DZ in the isogradient technique. Conclusive in this is 
that much less environmental data need be collected if a 
isogradient program were used. 
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VI. COMMENTS -AND -SSvSMMENDATIONS 

This thesis treated the NUTRACK III data as a comparator 
for analysis for both modes of STUTRACK I, but no conclusive 
evidence has been gathered to determine which is more 
accurate. The Range's Autotape system was to have determined 
this, but was so imprecise that a decision was unreachable 
without detailed analysis of the Autotape results. 

In all Z calculations at both stations the source depth 
plus Z was less than the array depth by three-to-ten feet. 

No hypothesis has been reached for the cause of the 
disappearance of these few feet. 

Further study of the data is needed to determine if the 
results from the intermediate stations are consistent with 
the conclusions reached here. Then, these conclusions should 
be tested using sound- velocity profiles representing other 
Dabob Bay oceanographic conditions. Special consideration 
should be given to the apparent inaccuracy in Z. 
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FIGURE 3. HYDROPHONE ARRAY 
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->l [<- 

RAWZ SIN(XTILT) 



The SIN(XTILT) is positive when the X transducer is 
above the horizontal plane passing through the C 
transducer, and negative when the X transducer is 
below the horizontal plane passing through C, 



FIGURE iJ. TILT CORRECTION 




FIGURE 5. RAY PATH REFRACTION 
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FIGURE 6. ISOGRADIENT THEORY 
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FIGURE 7. VELOCITY PROFILE 
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FIGURE 8. VELOCITY PROFILE 
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FIGURE 9. TRANSMISSION LOSS 
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FIGURE 10. RBNX, Station 1-25 
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FIGURE 11. RBNY, Station 1-25 
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FIGURE 12. RBNZ, Station 1-25 
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FIGURE ; 13. SLANT RANGE, Station 1-25 
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FIGURE 14. RBNX, DRIFT CORRECTED, Station 1-25 
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FIGURE 15. RBNY, DRIFT CORRECTED, Station 1-25 
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FIGURE 16. SLANT RANGE, DRIFT CORRECTED, Station 1-25 
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FIGURE 17. RBNX, Station 6-50 
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FIGURE 18 . RBNY, Station 6-50 
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FIGURE 19. RBNZ, Station 6-50 
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FIGURE 20. SLANT RANGE, Station 6-50 
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FIGURE 21. SLANT RANGE, DRIFT CORRECTED, Station 6-50 
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FIGURE 22. PHI, Station 1-25 
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FIGURE 23 . THETA, Station 1-25 
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FIGURE 24. PHI vs. THETA, Station 1-25 
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FIGURE 25 . PHI, DRIFT CORRECTED, Station 1-25 
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FIGURE 26.. .THETA, DRIFT CORRECTED, • Station 1-25 
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FIGURE 27. PHI vs. THETA, DRIFT CORRECTED, Station 1-25 
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FIGURE 28. PHI, DRIFT CORRECTED, FROM COR COORDINATES, Station 1-25 
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FIGURE 29 . THETA, DRIFT CORRECTED, FROM COR COORDINATES 
Station 1-25 
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FIGURE 30, PHI, Station 6-50 
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FIGURE 31, THETA, Station 6-50 
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FIGURE 32. PHI vs. THETA, Station 6-50 




n n j 

till 



71 



FIGURE 33- DX, Station 1-25 
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FIGURE 34. • DY, Station 1-25 
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FIGURE 35. DZ, Station 1-25 



ft* 




7 J l 



FIGURE 36 . DX, Station 6-50 



FT'T' 




75 



FIGURE 37. DY, Station 6-50 
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FIGURE 38. DZ, Station 6-50 
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APPENDIX A 



STUTRACK I COMPUTER PROGRAM 
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CALL FIX(SCX) 

CALL HI STF(DXIyNfNBAR) 
CALL FIX(SCY) 

CALL HI S T F ( DYI, N,NBAR) 
CALL FIX(SCZ) 
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